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Abstract 

We study the appearance of bound states in the spin gap of spin- 1/2 ladders induced by weak bond disorder. Starting 
from the strong-coupling limit, i.e., the limit of weakly coupled dimers, we perform a projection on the single-triplet 
subspace and derive the position of bound states for the single impurity problem of one modified coupling as well 
as for small impurity clusters. The case of a finite concentration of impurities is treated with the coherent-potential 
approximation (CPA) in the strong-coupling limit and compared with numerical results. Further, we analyze the 
details in the structure of the density of states and relate their origin to the infiuence of impurity clusters. 
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During the past decade spin ladder systems have 
attracted large interest among theoretical and experi- 
mental physicists (see Ref. [1] for a review) leading to 
considerable progress in understanding properties of 
such systems. In this paper we focus on spin ladders 
with weak bond-disorder and we allow couplings along 
both rungs and legs to be modified. The Hamiltonian is 

JV 

H = Y.^Jk,iSi,^-Si,2 + {Jl,i.iSi.i-Si+i,i + 1^2)] .(1) 
1=1 

A'^ is the number of sites and 5*;,^ denotes a spin- 1/2 
operator acting on site I on leg i; i — 1, 2. As we have 
shown in Ref. [2] the presence of bond-impurities can 
induce bound-states in the spin gap. In extension of 
our previous work [2] we address two examples and 
present results for the density of states (DOS) of the 
single-triplet excitation: (i) a random, binary distribu- 
tion of impurities on rung sites with a large impurity 
concentration and (ii) binary distributions of modified 
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couplings on legs and rungs. Thus, impurities are dis- 
tributed according to the probability density 

P{Ji) = cr,l5(J, - Jr,l) + (1 - cb.,l)S(Ji - Jr,l) (2) 

where Jb.,l denote the couplings of the pure system 
and Jr L are the modified couplings. In both cases, we 
assume the perturbations to be uncorrelated. Further- 
more, the concentrations for rung impurities cr and for 
leg impurities cl are independent. 

To treat this problem, we proceed by first mapping 
the spin-1/2 operators onto bond-boson operators sj 
and J [3], where sj creates a singlet on the Ith. rung 

and ; a triplet excitation with orientation a = x,y, z, 
respectively. Next, the singlet is integrated out and we 
restrict our study to the one-triplet subspace. Thus, 
triplet interactions and quantum fiuctuations are ne- 
glected which is equivalent to first-order perturbation 
theory in Jl / Jr for the pure model. 

The resulting effective Hamiltonian reads 

//off = El-^R.'i^' + ^^^"^"^^^(tl+.ti + H.c.)] ; (3) 
;=i 

omitting the spin index a but keeping in mind that 
the spectrum is threefold degenerate. The dispersion 
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of the homogeneous model is ek — Jr + Jl cos{k) , k 
being the momentum (see, e.g., Ref. [4]). The spec- 
trum of this effective model in the presences of impu- 
rities can be obtained numerically for very large sys- 
tem sizes and has to be sampled over many realizations 
of impurity-distributions (numerical-impurity averag- 
ing, NAV). Analytically, the self-energy T.{E) (E 
E + iO"*") of the one-triplet Green's function G{E) is 
accessible via diagrammatic techniques. Here we ap- 
ply the coherent-potential approximation (CPA) to the 
case of cr > 0; cl = 0; the self-energy being a self- 
consistent solution of the equation (see, e.g., Ref. [5]) 

Cr (5Jr 



E(£) = 



(4) 



1 - G{E) [dJR - S(S)] • 

(5Jr is defined as Jr — Jr, compare Eq. (3). 

Figure 1 contains our results for the DOS. First, 
we discuss the comparison of our numerical data and 
analytical computations for case (i), i.e., a binary dis- 
tribution of rung-impurities with Jr = 1, Jl = 0.2, 
and (5Jr = —0.2 as shown in Fig. 1 (a). The concentra- 
tion of impurities is cr — 0.5. Note that leg couplings 
are not perturbed. Both curves (solid: NAV, dashed: 
CPA) are symmetric with respect to i? = 0.9 while the 
original band for cr = was centered around E = 1. 
In the limit of cr = 1 the band will be centered around 
E — 0.8 possessing a band width of 0.2. The CPA gives 
a qualitatively good description of the overall struc- 
ture of the DOS. As the effective Hamiltonian (3) is 
that of a one- dimensional non-interacting disordered 
system, one expects a smooth DOS in the thermo- 
dynamic limit. However, the NAV shows additional 
structures [see, e.g., the peaks a,. . . ,h in Fig. 1 (b)] 
which are stable against variation of N for the system 
sizes investigated. These peaks can be attributed to 
the presence of impurity clusters (see below). Small 
impurity clusters like, e.g., (RR) (two modified cou- 
plings on neighboring rungs) or (ROR) (two modified 
couplings on rungs separated by one non-perturbed 
rung) appear as soon as the concentration is finite and 
lead to discrete peaks in the DOS. The corresponding 
cigcu-cncrgios can be computed exactly by solving 
Schrodinger's equation (see Ref. [2]). 

The influence of small impurity clusters will be fur- 
ther elucidated in the discussion of the second example 
(ii) where we have rmmerically computed spectra in 
the presence of fiintc concentrations of both rung and 
leg impurities. Fig. 1 (b) shows the resulting DOS for 
Jr = 1, Jl = 0.2, 5Jr = -0.2, 5Jl = J[ - Jl = 0.3 
and Cr — 0.1, Cl = 0.05. Note that cl = 1 corresponds 
to all leg couplings on both legs being modified. 

In Fig. 1 (b), the shape of the original band is still 
visible for 0.8 < E < 1.2 as the impurity concentrar- 
tions are fairly small. Both below and above the band, 
bound and anti-bound states appear. Their positions 
can be related to certain impurity clusters in the fol- 
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Fig. 1. Panel (a): comparison of CPA and NAV for a large 
concentration of rung impurities (cr — 0.5, cl — 0). Panel 
(b): NAV for bond-impurities on both rungs and legs. The 
letters label peaks originating from certain impurity clusters 
as indicated in the legend. See text for more details. In both 
cases, NAV was performed over typically 5000 random samples. 

lowing way. We consider a single impurity or a single 
cluster in an otherwise clean system and compute the 
eigen-energies of the corresponding states outside the 
original band, if such states are present (see Ref. [2] for 
a detailed discussion). Then, the main peaks visible in 
Fig. 1 (b) can easily be identified as being caused by 
small clusters or single impurities such as (R), (RR), 
or (RL) . Here, (RL) denotes a modified rung coupling 
and one of the neighboring leg-coupling modified. 

Thus, at low concentrations, the structure of the 
DOS can basically be explained in terms of non- 
interfering, isolated impurity clusters. The physical 
reason is that the wave-functions of bound-states are 
the less extended in real space the larger the ratio 
SJk/Jr, or SJl/Jl, respectively, is (see Ref. [2] for 
more details on the single- impurity wave- function). 
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